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Abstract. Phase transitions arc ubiquitous, exist in all fields of science in 
one form or another. The most common example in condensed matter physics 
is the Landau thermal phase transition, which is from a liquid to a solid phase 
or vice versa. Here, we systematically explore and develop unequivocal the- 
oretical strategies, going beyond the total-energy minimization techniques to 
understand what constitutes the thermal phase transition. We prove the ex- 
istence of finite-temperature continuous quantum phase transitions (^QPT) 
during solidification and melting processes such that ^QPT is responsible for 
all first-order thermal phase transitions. In fact, this ^QPT is related to chem- 
ical reactions where the quantum fluctuations are caused by thermal energies, 
and it may occur maximally for temperatures much higher than zero-Kelvin. 
To extract the quantitative information related to ^QPT, we use the ion- 
ization energy theory to derive the electronic-excitation and atomic-disorder 
(or symmetry-breaking) entropies. Subsequently, we exploit the energy level 
spacing renormalization group method to renormalize (i) the Bose-Einstein 
distribution and (ii) the specific heat capacity. We also conclude that the 
proofs developed herein can lead us to the unification of thermodynamics and 
quantum mechanics such that one can evaluate the changes in specific heat 
due to wave function transformations (quantum mechanical effects) when one 
crosses over from a solid to a liquid or vice versa at a constant temperature. 



PACS: 05.30.Rt; 05.30.Fk; 05.70.Jk; 05.70.Ln 
MSG: 82C26; 82C27; 82C22 

1. Introduction 

Unlike any thermal fluctuation, the fluctuation at a critical point of a flnite- 
temperature continuous quantum phase transition (^QPT) is specifically called 
quantum fluctuation (even for flnite temperatures) because it is related to the 
properties of electrons at a constant temperature, and therefore, it is still subject 
to the established Heisenberg uncertainty principle. An example in this respect 
is the existence of a ^QPT (due to changing pressure), which is well-known for 
ice at a constant temperature (0 < T < 273. 16K). The readers are referred to an 
elementary and straightforward discussion on these two issues (quantum fluctuation 
and jJQPT in ice) given by Sachdev [T]. But little did we know then, that physically, 
and logically, not all quantum fluctuations (wi^2, - ) will have to become "irrelevant" 
when T > OK, i.e., not all uj — > maximum when T — >■ OK. If the statement — 
all quantum fluctuations become maximum when T OK is true, then chemical 
reactions cannot exist. The above argument can be proven if we can prove the 
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existence of some quantum fluctuations tliat become large for finite temperatures 
or when T > OK. 

The quantum fluctuations discussed in Refs. [5J |3] are due to temperature- 
independent external tuners, namely, magnetic fields, pressure and chemical compo- 
sition affecting the spin carried by the electrons. Therefore, the quantum fluctuation 
discussed therein are entirely related to spin ordering and spin disordering, giv- 
ing rise to antiferromagnetism (AFM), ferromagnetism (FM) and paramagnetism 
(PM). Even if the system is none of these magnetism-related, spin will still come to 
play its role in one way or another when T — >■ OK. For example, a Cooper pair [3] 
in a superconductor consists of two electrons, one with (spin-down) and the 
other with (spin- up). Anyway, it is well known that the strength of this spin 
ordering (AFM or FM) or Cooper pairing (pairing of electrons with spin up and 
down) is maximum for T — OK, and the ordering or pairing gets destroyed when 
the electrons got excited, or when the Cooper pairs break up for some T > OK. 

Rightly so, Sachdev has correctly pointed out that the quantum phase transitions 
discussed in Refs. [Ij i2j do not include the temperature effect. Furthermore, these 
quantum fluctuations can be examined even for T > OK provided that uj/T can be 
expanded for higher temperatures if hu > k^T. Therefore, we can certainly add 
here that the temperature effect was excluded (or included as a form of expansion on 
uj/T) because the objective there and then was to evaluate the qualitative changes to 
the ground-state wave function of a particular system by varying the temperature- 
independent external tuners (for T = OK), not because QPT has been proven not 
to exist for T ^ OK. As a matter of fact, we prove here the existence of a QPT for 
r> OK. 

Now, to prove the existence of a particular QPT for T ^ OK, we need to un- 
derstand how the thermal energy initiates a quantum fluctuation during melting 
and solidification processes, and then evaluate the relationship between the ther- 
mally driven first-order TPT and ^QPT. Here, ^QPT requires quantitative changes 
to the properties of the electrons (due to the changes in their energy levels) that 
need some qualitative and/or quantitative changes to their wave functions. These 
changes are driven by thermal energy, not by the temperature-independent external 
tuners listed above. In contrast, the commonly accepted quantum phase transition 
for T = OK (okQPT) also requires quantitative changes to the properties of the 
electrons (due to the changes in their energy levels) that need some qualitative 
and/or quantitative changes to their wave functions. However, the above changes 
for T = OK are driven by temperature-independent external tuners alone, not by 
thermal energy. Hence, we should not reserve all "quantum phase transitions" 
exclusively for transitions driven by temperature-independent external tuners for 
T = OK, and not for transitions driven by thermal energy. 

Such a reservation is somewhat naive because the external drivers can be any- 
thing, including thermal energy, for as long as these drivers can initiate quantum 
fiuctuations for any T > OK. We emphasize here that a particular quantum fiuc- 
tuation that causes a phase transition is the one that should determine whether a 
given phase transition has a quantum origin or not, not by some ad hoc arguments. 
Of course, any given QPT can be driven by thermal energy, with further assistance 
from the temperature-independent tuners, or a QPT can be suppressed if one ex- 
ternal tuner competes with other tuners. In this case, it is just a matter of finding 
which tuner is in command of which interaction and quantum fiuctuation. The 
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quantum fluctuation we are concerned with here is the fluctuation in the energy- 
level spectra (energy-level spacing) at a critical point. These energy levels refer to 
electrons, of course. In any case, we have said enough. At the end of this section, 
some of you may come to realize that the quantum-fluctuation causing thermal 
energy holds the key to unify thermodynamics and quantum mechanics. 

Having explained that, we now continue with our introduction to ^QPT. The 
quantum fluctuation associated to ^QPT occurring at a flnite temperature quantum 
critical point (tQCP) requires quantitative changes to the wave functions, such that 
the quantum fluctuation with an energy scale, huj can be written as, 



where h is the Planck constant divided by 27r, fee denotes the Boltzmann constant, 
uj{y) and Tg are the doping (y)-dependent frequency and some critical temperature 
(a rQCP), respectively. Later, we will come to see why the wave functions need 
to be changed quantitatively for ^QPT. For a okQCP, the relevant inequality is 
simply hjjj(]ji) < huj{y2) < huj{yy,) because Tg — OK. Here, j/2 and j/3 are related 
to different chemical composition, while 6 represents a certain physical property 
under investigation. 

For example, for the melting process (chemical dissociation), Tg = Tp^^^^j'"®, and 
it satisfies the classical phase transition (due to thermal fluctuation with increasing 
temperature for a given chemical composition), without any microscopic under- 
standing as to the reasons why and how the translational symmetry breaks (ice to 
water). This does not imply melting is a classical process not involving electrons. 
For instance, we have proven how the energy level spacing evolve and give rise to 
a changing interaction strength between a polarizable water molecule physisorbed 
onto the MgO surface, which then initiate the water molecule dissociation at a con- 
stant temperature [S] . Within the ionization energy theory (lET) [HI El [H] ; we can 
write hjjj{y) = S.{y) where ^(y) is the ionization energy of a particular compound, 
and it is also known as the energy level spacing within the renormalization group 
method [3]. Hence, during chemical dissociation with changing chemical composi- 
tion, ^{y) is a dynamical parameter, giving rise to quantum fluctuation. 

Importantly, much below the melting point, the fluctuation in ^(y) is caused 
by the energy level crossings, changes in the valence states of the constituent el- 
ements, and defects due to different chemical compositions (doping) along a con- 
stant temperature line. Such fluctuations (due to defects and valence states) are 
responsible for the changes in the constant-temperature and composition-dependent 
(i) resistivity of ferromagnets [10], (ii) resistivity in the normal-state of super- 
conductors [HI im HI] , (iii) non- blinking effect in quantum dots [13] and (iv) re- 
versible carrier- type transition in oxides [14] . Apart from that, experimental and 
theoretical developments of the popular zero temperature quantum phase transi- 
tions (okQPT) in strongly interacting solids were also investigated independently, 
namely, in two dimensional systems |15|, I16j , in solids with respect to spin order- 
ing [HlilSKinKinillTllSll^lMlllSlllS], in Bose-Hubbard model [27], and in 
antiferromagnets [551 [2H]- 



(1.1) 
(1.2) 
(1.3) 



fi^iyi) < ksTg, 
Sa;(y2) = k^Tg, 
^(2/3) > ksTg, 
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In the subsequent section, we introduce and explain why and how the molten 
alkali halides give rise to the first-order thermal phase transition, and how this tran- 
sition is related to the finite temperature quantum phase transition. Subsequently, 
we venture deeper into the details of quantum phase transitions and quantum fluc- 
tuations with respect to the ionization energy theory, which involves wave function 
transformations and electronic phase transitions. Along the way, the theory is 
discussed and crosschecked whenever necessary against the experimental observa- 
tions. Finally, we will conclude that (a) the finite-temperature continuous quantum 
phase transition is responsible for all first-order thermal phase transitions, and (b) 
at a critical point of a first-order thermal phase transition, a non-divergent time- 
dependent specific heat can be proven to exist. Both points in (a) and (b) involves 
some physics beyond the traditional phase transitions and critical phenomena. 

2. Further analyses I 

Here, we will give unambiguous logical proofs (supported by experimental results, 
of course) on the origin of thermal and quantum fluctuations due to temperature 
and doping (or £,{y)) within lET. But before doing so, we need to justify the reasons 
why we have selected salts made up of single-valent cations and anions from group 
17 of the periodic table. 

2.1. Thermal and quantum fluctuations. Table [T] lists the diatomic bonding 
energies and melting points of these salts such that the melting points decrease 
systematically with decreasing bonding energies from F to I. Such a decreasing 
trend can be precisely overlapped with the decreasing ionization energies of the 
anions, namely, Cf+(1681 kJmor^) > Cci+(1251 kJmor^) > CBr+(1140 kJmol-i) 

> ^j+(1008 kJmol"^). The values in the inequalities are given in Table [2] follow the 
values marked with " * " . Interestingly however, for a given anion, the decreasing 
melting points do not agree with the ionization energy values of the cations: Li, 
Na, K and Rb (follow their values listed in Table [2] marked with " f "). 

For example, the mehing points for NaF (996°C) > KF (858°C) > LiF (848°C) > 
RbF (833°C) cannot be overlapped with ^u+{520 klmor^) > ^Na+(496 kJmol-i) 

> ^K+(419 kJmol"^) > CRb+(403 kJmol^^). This is not surprising because melt- 
ing is a process directly proportional to bonding strengths (see Table [T]), which 
have been discussed earlier [T3] within lET where large ionization energy values 
of cation-like ions (^0=+ < ^0^+) in a given molecule (C^+0^~ or 0^j^0^2)) do 
not necessarily imply stronger bonds. Fortunately, we have proven that a stronger 
bond is predictable from lET if one considers the ionization energies of anions (say 
O in H2O), because oxygen defines the ability to attract electrons from an atomic 
hydrogen [S]. In this case, smaller means weaker O— H bond, which is consistent 
with the above-stated ionization-energy inequalities for anions (group 17 elements) 
and the melting points. 

In Table [2] we have deliberately selected systems consisting of single valent 
cations and anions from group 17 to avoid effects from different electronic interac- 
tions due to different number (and type) of constituent atoms in a given molecule. 
To understand this point, we again use molecular systems, in which, for a NO2 
molecule, N acts as a cation, while O as an anion, which means that we need to 
consider the valence state of 4-1- since four electrons have been transfered from N''+ 
to Oj". This electron-transfer is due to ^^4+ < ^0*+ (see Table [l]). In contrast, N 
is an anion in NH3 molecule due to > where three electrons from three 
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Table 1 . Experimental values for the melting points and diatomic 
bonding energies of salts obtained from Refs. [30| and [31], respec- 
tively. The systematic decrease of the melting points and diatomic 
bonding energies with respect to anions (from F to I with increasing 
Z) satisfy the decreasing ionization energies for the same anions, 
from F to I. See the text and Table [2] (follow the values marked 
with " * ") for details. 



Salts 


IVlGlting points 


Diatomic bonding 




(°C) 


energies (kJniol^^) 


LiF 


848 


577 


LiCl 


610 


469 


LiBr 


552 


419 


Lil 


469 


345 


NaF 


996 


519 


NaCl 


801 


412 


NaBr 


747 


367 


Nal 


660 


304 


KF 


858 


498 


KCl 


771 


433 


KBr 


734 


380 


KI 


681 


325 


RbF 


833 


494 


RbCl 


715 


428 


RbBr 


682 


381 


Rbl 


642 


319 



hydrogen atoms are transfered to nitrogen, giving N'^~Hj^. Thus far, the analyses 
are correct. However, the polarizability of the molecule NO2 does not solely depend 
on these four electrons contributed by nitrogen if we compare N^'^+Oj" with N'^~Hj^ 
because we cannot simplify the analysis by comparing ^n*+ with (3 electrons 
contributed by three hydrogen ions) only. If wc do so, then obviously we have 
Cn*+ > '?H+ that falsely allows us to conclude a^'~'^ < a^^^ because oc exp(— ^) 
where ad is the displacement polarizability. 

In other words, we cannot use the above inequality (^) to directly compare ad be- 
tween NO2 and NH3 molecules because N acts as a cation in a NO2 molecule, while 
it is an anion for the molecule, NH3. This implies that we need to take both the 
cationic- and anionic-effect into account explicitly for an accurate logical analysis, 
which have been correctly invoked in Ref. [S]. However, the anionic-effect can be 
neglected when we compare CO with O2 because in these molecules, oxygen is the 
anion and therefore a^'-' < is valid [14]. On the other hand, the justification 
required to neglect the cationic-effect for doped-Pnictide superconductors is given 
in Ref. [H]. Finally, taking both cationic- and anionic-effect into account means 
that for large (many-electron) cations such as K (Z = 19) and Rb {Z = 37), there 
will always be some significant amount of polarization (due to large screening) from 
the outer core electrons, even though effectively, and Rb"*" are single-valent ions. 
This second- or third-electron polarization is negligible for few-electron atoms (due 
to small screening), namely, Li {Z — 3) and interestingly, also for Na {Z — 11). 
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Table 2. Averaged atomic ionization energies (^) for individual 
ions and their respective valence states ordered with increasing 
atomic number Z. All the experimental ionization energy values 
were obtained from Rcf. 31 . 



Elements Atomic numbers Valence ^ 





Z 


states 


(kJmol~i) 


H 


1 


1+ 


1312 


Li 


3 


1+ 


520t 


N 


7 


1+ 


1402 


N 


7 


4+ 


4078 





8 


1+ 


1314 





8 


2+ 


2351 





8 


4+ 


4368 


F 


9 


1+ 


1681* 


Na 


11 


1+ 


496t 


CI 


17 


1+ 


1251* 


K 


19 


1+ 


419t 


Br 


35 


1+ 


1140* 


Rb 


37 


1+ 


403t 


Sr 


38 


2+ 


807 


I 


53 


1+ 


1008* 


La 


57 


3+ 


1152 



Note here that we need smaller and larger ionization energies for cations and an- 
ions, respectively for stronger ionic bonds. Therefore, incorporating this additional 
amount of polarization (for cations) implies larger ^cationic or equivalently, smaller 
^anionic that should reduce the melting point as has been observed experimentally 
(see Fig. [T]) . Importantly, this correction does not apply for anions if a given system 
consists of mostly ionic bonds because the anions are judged solely on their ability 
to attract electrons from the cations. On the contrary, for systems with mostly co- 
valent bonds, we need to consider both cationic- and anionic-effect simultaneously 
as carried out in Ref. [5]. The existence of this additional polarization is consis- 
tent with the experimental results shown in Fig. [T] when one compares the slope, 
dTfi/d^ for salts containing Li and Na, for instance (dTg/d^jLi ~ (dr6(/d^)Na, while 
(dTg/d^K « (dre/d.^)Rb (follow the solid lines in Fig.[l]). 

2.2. Finite-temperature continuous quantum phase transition in a two- 
atom system with asymmetric polarization. To see why and how the above 
quantum phase transitions (due to changing T and chemical composition) come to 
play to give rise to two different phases (solid and liquid), a model Hamiltonian 
would be helpful for some readers, though it is not necessary. Unfortunately, we 
cannot write a quantum mechanical Hamiltonian right at the melting point because 
both the Hamilton operator and the wave functions are time-dependent, and they 
change quantitatively with time [35] . Instead, our strategy here is to first evaluate 
a Hamiltonian consisting of two identical atoms with interactions. After obtaining 
the attractive interaction energy between these two neutral polarizable atoms, we 
will then add the polarization induced attractive interaction term appropriately to 
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Figure 1. Melting points of Li(I,Br,Cl,F), Na(I,Br,Cl,F), 
K(I,Br,Cl,F) and Rb(I,Br,Cl,F) salts are plotted against the an- 
ionic (I, Br, CI, F) ionization energies, which indicate a direct pro- 
portionality between them (follow the solid lines). INSET: The 
vertical down-arrow (solidification) denotes a first-order thermal 
phase transition (from liquid to solid) for LiCl when the tem- 
perature is reduced from T > 610°C to T = 552°C. The hori- 
zontal arrow pointing left captures the continuous (purely elec- 
tronic) quantum phase transition (from solid to liquid) at constant 
T ~ 552°C by changing the chemical composition systematically, 
LiCli_aiBrai-a2la2 such that LiCl — > LiBr — > Lil, = 1 — oi, 
2/2 = oi - a2 and yi — 02. 



Eq. |2.4| to obtain the total attractive interaction between two atoms with asym- 
metric polarizability. This total attraction is shown to exist just above the melting 
points of the single- valent ionic bonded salts listed in[l] The first step is to consider 
two single-electron, identical and neutral atoms (A and B) depicted in Fig. [2F A) . 
The relevant Hamiltonian (with polarizable electrons, ea and ee) is given bv |33j 



(2.1) 



H 

Ho = 
Hi = 



Ho + Hi 



2m 
1 



krl 
2 

MM 

R 



2m 2 ' 



i-e){e) 

R-TA 



(-e)(-e) 
i? - TA + re 



which can be separated into two harmonic oscillator Hamiltonian after taking R ^ 
''A,B where ta and re are the coordinates for electrons, ca and cb, respectively, [33], 



(2.2) 



H = 



1 

2 m 



P4 



k - 



2neQR^ 



1 
2m 



k 



e 

27reoi?3 



To obtain Eq. 2.2 we have carried out the substitutions listed below 
1 



(2.3) 



1 



^^(^A ± 7-b) and p± = -^{pa±Pb)- 
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Figure 2. (A): In writing Eq. 2.1 and its solution in Eq. 2.4 



we have actually assumed the two sketched atoms (A and B) to 
be polarizable, neutral and identical with discreet energy levels. 
Moreover, each atom, A and B has one equally-polarizable elec- 
tron, namely, ca and cb, respectively. Here, the negatively charged 
electrons and positively charged nuclei are arranged in this man- 
ner, + — + — along the same axis. (B): On the other hand, 
^Coulomb ^^"^ rcnormalized screened Coulomb potential between 
atomic A and B, which has been added by hand into Eq. 2.5 (the 
first term on the right-hand side). This particular term (V'(^(7u°omb) 
imposes the condition that both A and B are also single electron 
atoms, but with a little twist such that atomic A has the least po- 
larizable electron (ca), while atomic B has the easily polarizable 
electron (cb). In the text, we will also denote atomic A as the an- 
ions (namely, atomic CI), whereas atomic B represents the cations, 
namely, atomic Li. The charge carried by each nucleus, A and B 
in both sketches, (A) and (B) is +e. 



Subsequently, one obtains the ground state energy for the above Hamiltonian |33j 



(2.4) E = K\ k-{eV2^e,R^) ^ k + (^j2^^ 

2 V m V m 

Here, pk,b are the electrons momenta, m is the electron mass, h is the Planck con- 
stant divided by 27r, R is the distance between two neutral and identical atoms and 
eo is the permittivity of free space. The bounded electrons are labeled as ca and Cb , 
which are bounded within nuclei A and B, respectively. In this simple Hamilton- 
ian, the electrons (ca and Bb), and the nuclei (A and B) carry the charge, — e and 
-|-e, respectively. Furthermore, we have denoted Hq as the ordinary Hamiltonian 
that takes the kinetic energy terms, and two semiclassical harmonic-oscillator type 
electron-ion potential terms into account (see Eq. |2.1[ ). On the other hand, we also 
have the interaction Hamiltonian denoted by Hi to capture all the Coulomb inter- 
action between atomic A and B (see Eq. 2.1 and Fig. [2f^A)). We stress here that 



neither Hq nor Hi takes the induced repulsive Coulomb forces due to polarized 
electrons into account. However, the attractive Coulomb interactions between the 
valence electrons and their respective nuclei have been correctly included in Hq. 

In our previous work ^7] , we have proven the existence of a blue-shifting van der 
Waals (vdW) attractive interaction between two neutral atoms, after carrying out 
the energy-level spacing renormalization on the interaction potential constant, fc, 
and after imposing an additional condition that requires atomic A to be least polar- 
izable, while atomic B as the easily polarizable one (see Fig.[2|^B)). This condition 
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gives the asymmetric polarization that leads to a stronger attractive interaction be- 
tween atoms. Such a stronger attractive interaction can lead to chemical reactions 
as proven in Ref. The stronger attractive vdW interaction energy [37] 



(2.5) 

(2.6) E- 
and from Ref. C 



Coulomb 



E 





(2.7) 

From [5] 
(2.8) 



T T-c— ion 
Coulomb 



47reo|RA - tb 













j jexp 




+ exp 




} 



k — k exp [A^] , 



where E and k are the renormalized energy (from Eq. (2.4)) and interaction poten- 
tial (spring) constant, respectively, fi is the screening constant of proportionality, 
A = (127reo/e^)as, a_B is the Bohr radius of atomic hydrogen, Ra and re are the 
coordinates for the nucleus A and electron ee, respectively. Here, h{k/mY/'^ = tooj 
too is the renormaliz ed e nergy when Hi = 0, a nd V^^^i^{S,) is negative guaranteed 
by l/\/2— 1 in Eq. (2.6) and — e in Eq. (2.7). The term on the right-hand side 
of Eq. (2.6) assumes both atoms, A and B are identical, hence their electrons are 



equally polarizable. In Eq. (2.7 1, we have a unidirectional electron-ion attractive 
interaction between the easily polarizable ee and the nucleus A, screened by ca via 
^A (see the term in the curly bracket in Eq. ( |2.7| )). The imposed asymmetric polar- 
izability is captured by V^~uYomb that the respective atomic A and B represent 
the anions (least polarizable) and cations (easily polarizable), which is required to 
understand the transition from a liquid to a solid phase. 

To understand the types of phase transition in the above model Hamiltonian, one 
needs to observe the evolution of the Hamiltonian from Hq (non-interacting gas) 
to Ho -\- Hi (weakly interacting liquid) to Hq + Hi-\- V^couiomb (strongly interacting 
liquid) to ifgoiid (solid). Therefore, one can actually invoke two types of phase 
transitions from Eqs. ( |2.1[ ) and (2.7), from a gas to a liquid phase, and then toasolid 
phase. For example, taking Hi 



and V^Q^Yo'inb ~^ 0, one immediately obtains the 
non-interacting gas phase. For a liquid consisting of identical atoms, and if the atom 
is least polarizable, then their only interaction is the long-range vdW attraction such 



that Eq. (2.1) is a sufficient model Hamiltonian for this weakly interacting liquid. 
On the contrary, if the liquid consists of two different atoms such that one of them 
is highly polarizable (cation- forming (or electron-donating) atoms) , while the other 
atom contains highly non-polarizable electrons (anion-forming (electron-accepting) 
atoms) , then the liquid consisting of these atoms forms a strongly interacting liquid. 
For example, any molten solid close to the solidification point, or slightly above the 
melting point. 

Now, you should not be carried away into thinking that a Coulombic liquid is 
different from the atomic molten system. For example, a system that consists of 
highly polarized atoms forms a strongly interacting Coulombic liquid, while the 
weakly interacting Coulombic liquid consists of weakly polarized atoms. Obviously, 
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unpolarized or uncharged system cannot form such a (Coulombic) hquid. Secondly, 
the model discussed above nicely corresponds to the experimental results listed in 
Tables [T] and [2j namely, LiF, NaF, KF and RbF have the largest melting points 
relative to other listed anions because the ionization energy for F is the largest (with 
the strongest ability to attract electrons), compared to I, Br and CI. This implies 
that at the melting point, the breaking and forming of stronger bonds will lead 
to a strongly interacting Coulombic liquid, as it should be. Finally, based on the 
above model, the bond-formation and-breaking mechanism is due to the quantum 
phase transition (giving rise to electronic phase transition), which requires a wave 
function transformation. The electronic phase transition stated above has got to 
be nonadiabatic because of electron transfer during bond-formation or -breaking 
(chemical reactions). However, the adiabatic approximation can still be applied, if 
and only if the Hamiltonian (-ff (i)reaction') describing this chemical process is known 
explicitly [35]. 

It is well-known that when one cools a molten solid, then one needs a specific 
amount of thermal energy or simply heat (Qmciting) to initiate the quantum fluctu- 
ation right at a critical point (at the melting point) to give rise to a wave-function 
transformation (electronic phase transition) 'S^J . Any heat larger than Qmciting will 
not activate the wave function transformation required for the liquid to solid phase 
transition. For example, any Q > Qmciting implies Q > {V^~^^°^^^) and therefore. 
Coulomb repulsion dominates. On the other hand, if Q < Qmoiting then one has 
Q < which implies that the heat is insufficient to overcome the bonding 

energy in the solid phase, prior to melting. In the remainder of this paper, our focus 
is to understand how Q is related to entropy, specific heat and to the energy-level 
spacing (^) of the molten system. Our strategy here is to prove that Qmciting drives 
the changes in i^, which in turn initiates the electronic phase transition. Within 
the Hermitian quantum mechanics, the said electronic phase transition is caused 
by the wave function transformation [321. Therefore, a fluctuation in ^ corresponds 
to a fluctuation in the many-body wave function, and subsequently, is responsible 
for QPT. In other words, any external tuner, including thermal energy can cause 
a fluctuation in ^, which implies a quantum fluctuation that eventually could give 
rise to a QPT. Later, we will expose why and how the fluctuation in ^ (at a critical 
point) forces us to impose time-dependence into ^. 

3. Renormalized entropy and Bose-Einstein distribution 

Figure [l] shows the direct proportionality between melting points and ionization 
energies of the constituent anions (I, Br, CI, F). The respective vertical and horizontal 
arrows in the inset of Fig. [l] indicate the first-order TFT and ^QPT. In order to 
expose the existence of these phase transitions, we need to recall the first and second 
laws of thermodynamics. These two laws can be combined to obtain 

(3.1) dU = SQ + 6W ^ TdS - PdV, 

in which, the change in the internal energy, dU of a closed system equals the 
amount of heat (Q) absorbed and the amount of work {W) done by that system. 
Here, S is not an exact differential because the changes in Q and W depend on 
the thermodynamic path (or independent of a particular system and process) , and 
therefore, Q and W are not state functions. Here, P and V denote pressure and 
volume, respectively and we consider only reversible processes. Unlike Q and 
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the thermodynamic variable U, P and V are state functions that are unique for a 
given system and process. 

The second law is given by, SQ = TdS where S is the entropy, another path- 
independent state function. All we need now is the relationship connecting 5 to ^ 
that can be obtained from the derivation of the ionization energy based Fermi-Dirac 
statistics (iFDS). Denoting N as the total number of particles with ni particles 
have energy (Eq ± n2 particles with energy (Eq ± and so on implies that 
N = ni + n2 + ■ ■ ■ + n,„. As a consequence, the number of ways for qi unoccupied 
quantum states to be arranged among ni particles is 



(3.2) r!(ni,gi) 



<Zi! 



ni!(gi - rii)! ' 



satisfying the condition that one empty quantum state can accommodate only one 
particle. Subsequently, the total number of ways for q quantum states {q = qi + 
• • • + Qm) to be arranged among N particles 



(3.3) r!(7V,g)=J] 



ni\{qi - ni)\ 



The most probable configuration for certain T can be obtained by maximizing the 
number of ways one can arrange particles in qi empty quantum states or, we 
need to maximize D.{N,q) subject to the restrictive conditions, 

(3.4) Y,n^ = N, ^dn, = 0, 



(3.5) ±O^^Eo±C^E, Y^iEo ± 0^dn^ = 0. 

i i 

These conditions impose that the total energy, E and the total number of particles, 
N for a given system are always constant. Here, Eq is the total energy for OK. The 
method of Lagrange multipliers |33j can be employed to maximize Eq. (3.3) by first 
introducing a new function, Fnew(' • ■ ; A*j '^b, •••) = •••+ fifi + Ab/2 + • • • , such 
that 

(3.6) F^,^ = \nn{N,q) + fiY^n, + XsY^i^n ^ 0^n^, 

i i 

where ni = n2 = ns • • • = = 1 to imply one empty quantum state (q) or energy 
level can accommodate only one particle (either electron or hole). Subsequently 
(we maximize -Fncw), 



Equation (3.7) can be solved subject to the conditions, 3> 1, S> 1, qi > rii 
that give us the license to use Stirling's approximation (In (rij!) w In n^ — rii), and 
after some algebraic rearrangements, one gets 

(3.8) 



q exp[^ + AB(i;o±0] +1' 
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which is the iFDS. Taking exp [fi + Xb{Eo ± ^)] » 1, ^ = -XbE^ and dividing all 
the terms by T will lead us to the energy-level spacing entropy 

(3.9) 5, = -(^°^^)+^g^-Vln^, 
^ ^ ^ T (AB,A)r g' 

in which, 

(3.10) Ab — for constant ^ and 

kbT 

(3.11) A = ^"^"^^ for constant T > 0. 

The derivation for Ab is well known in classical thermodynamics and is also given 
in Ref. [5] within the lET formalism, while the proof for A is available in Ref. [7] . 
Here, Ep denotes the Fermi level for T = OK, independent of any disturbance. 



The entropy S*^ given in Eq. (3.9) decreases logarithmically (S*^ -> — oo) when 
N/q — > (ordered) such that, the inequality < N/q < 1 strictly corresponds 
to — cxDijiin < < Oinax- Furthermore, within the set of real negative numbers, 



including zero (5^ £ M^) and Eq. (3.9), we have the correct correspondence between 
T -)■ and 5^ -> -oo (ordered). 

The other relevant entropy is due to atomic-disorder, Sj) — /cb InD where 

(3.12) D(^„, ,,) . n - n 

by taking S> 1, 3> 1 and defining No/qn as the ratio of the total number 
of identical atoms and empty lattice sites. Here, iVo — X^i^-i; 9D — X^i?*; 
A'b /go corresponds to the probability of excited identical particles from their lattice 
sites, which will tell us the magnitude of ionic conductivity. For example, if some 
large number of identical atoms can be excited, then this will lead us to large 
ionic conductivity, which is a telling sign that we are approaching the melting 
point of a particular solid. Anyway, one of the two restrictive conditions needed 
to maximize Eq. ( |3.12 ) remains the same (Eq. (3.4)) because the total number of 



particles (atoms) is also conserved. While the condition for the total energy has 



been renormalized (see Eq. (2.8)), 




(3.13) E[^^^'']."^-M/^=^ 

i 

to imply that the total energy (of the atoms, not electrons) is also conserved. 
Moreover, we have renormalized the total energy of the identical atoms (E), not 
of the electrons, to obtain E so that the electronic polarizability effect is taken 
into account. This means that the atoms are no longer considered to be the rigid 
vibrating ions, independent of their surrounding electrons. Somewhat similar to 
iFDS derived earlier, the new function, -F^^w this case is 

(3.14) =lnD(iVD,9D)+aEn, + ABE [£;e5^«]^ri„ 



where each n can take any value such that n £ N to imply one empty quantum 
state (q) or energy level can accommodate any number of particles. Here, N is the 
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set of natural numbers. Subsequently (we maximize F^^^), 

(3.15) ^ ^ ^[WM + . E ^ + Ab E [^e^^^ ^ = 0. 
^ ' dn, dn, ^ dn., ^ * 5n, 



We solve Eq. (3.15) by invoking ^ 1, ^ 1, qi > and the Stirling's approxi- 
mation (ln(ni!) sa n^lnni — rii). After some algebraic rearrangements, and taking 
a — Q for convenience because we are concerned with the changes to ^, we obtain 
the renormalized Bose-Einstein statistics (rBES) 

(3-16) ^ = , . 

* exp [\BEe2>^^] - 1 

The constant, a denotes the minimum energy prior to any excitation of atoms from 
their equilibrium lattice sites, somewhat similar to Fermi level for electrons. It is 
important to note here that go > -^d is required so that the probability is normalized 
to one, even if there can be any number of n particles allowed to occupy a single 
empty quantum state, g^. Similar to S*^, Su £ Mr . 

We are now ready to track the solidification of liquid LiCl from T2 > 610°C to 
Ti = 552°C as shown in the inset of Fig. [T] In this case, the dominant entropy 
change is due to the transition from the disordered (liquid state) to an ordered 
(solid) state in which, the liquid state is defined within 610°C < T < T2 and the 
solid state is bounded within Ti<T < 610°C. The value 610°C (Tg) is the critical 
point, or the melting point of LiCl and S'p"'""^ > 5*^°'"^ indicates the existence of a 
TPT, qualitatively that is. To claim liquid has a larger entropy than solid implies 
we have invoked an additional condition IS'^''"''^ - S^°^''^\ > ISp""'"^ - because 
the change of entropy due to electronic excitation is always smaller compared to an 
atomic-disorder induced entropy (due to broken translational symmetry) within a 
particular system. Of course, in the absence of this temperature-induced atomic- 
arrangement asymmetry, electronic disorder or is the dominant one. 

Apparently, the change in the entropy during the above-stated transition (solid- 
ification) is a first-order TPT because S^!^!^j.^ > S'D°T<Te such that there exist a 
nonequilibrium entropy (S'"™^ that contribute to this inequality, which is classically 
undefined at the critical point, T = Tg = 610°C where both sohd and liquid phases 
coexist. However, we can define S*"™ with respect to the time-dependent energy- 
level spacing, ^{t) and using Eq. ( |3.9| ) at the critical point due to time-dependent 
processes of breaking and forming of bonds at the solid] liquid interface (we 

will revisit this issue in detail shortly) for T = Tg. In this case, ^ continuously 
switches between ^''i^'d and for T — Tg. Consequently, one cannot write 

^sohd ^ 5*^°^^ < S'd''""^, which in turn implies that solidification and melting are 
first-order transitions. 

4. Renormalized specific heat 

Traditionally, a first-order TPT is defined to exist if there is a quantitative 
discontinuity in the thermodynamic variable. For example, the heat capacity (a 
thermodynamic variable) is discontinuous at the critical point due to discontinuity 
in the entropy itself as stated above. Here, the constant-volume heat capacity and 
entropy relationship can be obtained from 



dU 



dT 







^dT, 
T ' 
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using Eq. (3.1) and the definition, — dU/dT\v, in which Eq. (4.1) teUs us 
nothing about the changes in Cy during sohdification or at the critical point, except 
that it diverges }40j. Anyway, on the way to understand the origin of the first order 
TPT (given above), we have pointed out the existence of time-dependent ^(t), which 
fluctuates between ^'''i"''^ and ^"^'"^ for T = Te in the presence of both solid and 
liquid phases. Now, the above fluctuation in ^ can be associated to the existence 
of ^QPT at the critical point during the melting or solidification process in such 
a way that if we continue extracting heat from the LiCl system, then one has the 
liquid to sohd transformation due to [^"'i"''^ ^ f ""'i] > ^ C"^'"^] that freezes 

LiCl completely. The formation of bonds (or releasing of heat) during solidification 
occur continuously with respect to time that involves complete modification of the 
time-dependent many-body wave function (4'(i)). Within lET formalism, we do 
not require to know the changes in ^(t), instead, we have ^(t) as the fundamental 
interaction-strength functional, which changes continuously and originates from the 
chemical constituents of a system. In other words, the chemical composition at the 
solidjliquid interface changes with time, giving rise to this correspondence rule, 
^(i) ~^ ^{t)- For a given temperature (T 7^ Tg), if the chemical composition is 
i-independent (for example, in the solid or liquid phase) but r-dependent, then 
^(r) — > ^(r) where r is the electron coordinate. In our earlier work [Sj, we have 
renormalized Cy that can be exploited here to rewrite Eq. H.ll for T ~Tg 



(4.2) d5- (0 = ^!^dt = ^ exp 



which reinforces the logic that any change to the interaction strength during so- 
lidification, even though T-independent at the critical point, never ceases to be 
t-dependent, giving rise to S'uquid ~^ S'soUd and Ci'^"'^ ^ Ci^°"'^ transitions. These 
t-dependent transitions, S'^q^{t) and Cv{t)cqm can be readily captured via £_{t) 
without any divergence. In addition to that, we now claim that the solid ^ liquid 
transitions at T = Tg are associated to §^QPT. An essential point to note here is 
that we did not identify C„(T) in Eq. ( |4.2[ ) as the nonequilibrium specific heat at 
the critical point because Cy{T) reduces to C|,"i""^(T) during solidification, whereas, 
Cv{T) — ^ C^°'"^(T) during melting. We also remind the readers to take note here 
that unlike S'5 and S'd, {S,S} G M+ because {N/q^N^/qxi} e (0, 1]. To those who 
are not comfortable with {5*^, S'd} G Mr , we define an alternative equation 



(4.3) S^{+) = /CB In ^1 = In TV, 

exp + A(£'o ± C)J +1 



which can be obtained from Eq. (3.8) and recall here that N is the number of 
excited electrons. Similarly, 

in this case, is the number of ions being displaced (or excited) from their crys- 



by rewriting Eq. (3.8) such that 



tallographic sites. Apart from Eqs. (4.3 1 and (4.4 1, we can also enforce positivity 



(4.5) 
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because if Ni/q > N2/q and < iVi,2/g < 1, then 

(4.6) -^ln^>-^ln^, 

q q q q 

which will guarantee S'l > 0, 5*2 > and > S2 because N/q increases or decreases 
faster than \n.{N/q). For example, taking N — 1, q = x, f{x) = 1/x and g{x) — 
ln(l/x), one can write h[f,g] = (1/a;) In (1/a;), subsequently it is straightforward 
to show d^Z/da;"" > d^g/dx" is always true when x > 1 and m > 1 where x E M+ 
and m g N* . In addition, Si > S2 is consistent with increasing entropy if N/q gets 
larger, and consequently we will always have {S'j(+), 5^, 5d(+), 5,^} G where 
and N* are the set of real positive numbers and positive integers, including zero, 
respectively. 

In summary, any first-order TPT has got to go through a ^QPT at constant T 
such that 

(4.7) lim^QPT = TPT and ?QPT D TPT, 
h—>-o 

which means, any first-order TPT is a limiting case in the absence of quantum 
fiuctuations, and therefore j^QPT is a proper superset of TPT. These two claims 
are strong, thus we will need to provide logical proofs. The proofs consist of two 
parts that correspond to two claims. The first one is given above when we discussed 
the solidification phenomenon for LiCl, establishing the correctness of Statement 1. 

Statement 1; All first-order TPT must go through a thermal-assisted '^QPT 
at constant T. 



The second proof proves Eq. (4.7 1, which will be exposed shortly. 



5. FiNITE-TEMPERATURE CONTINUOUS QUANTUM PHASE TRANSITION 

Thus far, we have exposed the existence of ^^QPT during LiCl solidification at 
T — Tg (see the vertical arrow in the inset of Fig.nJ, which is commonl y acc epted as 



a first-order TPT without going into the details ofjiQPT. Proving Eq. (4.7) requires 
one to track j^QPT or the changes in £_{t) during solidification for constant T = Tg. 



Alternatively, one may also prove Eq. (4.7) by tracking the horizontal arrow pointing 



left shown in the inset of Fig.[T]by systematically changing the chemical composition 
via substitutional doping of CI with Br and followed by I such that LiCl — >■ LiBr — >■ 
Lil at constant T — 552°C, which is the melting point of LiBr. We will first address 
the ^QPT occurring during solidification (for constant T — 610°C). The above 
^QPT during solidification of LiCl system can be captured by isolating the system 
right at the critical point (Tg = 610°C) as shown schematically in Fig. |3] Here, the 
liquid phase Li and CI are indicated with filled red circles, while the same elements 
in the solid phase are drawn as blue circles, but we did not bother to completely fill 
the solid phase with blue circles. Relevant nonequilibrium processes occur at the 
boundary between these two phases. These processes are due to chemical reactions 
between the highly-polarized Li (small ^Li) and the least-polarized CI (large ^ci) 
giving rise to the chemical association between them, forming the solid phase, which 
releases energy as heat (— Q) into the liquid phase. This energy transfer increases 
the kinetic energy of the liquid-phase Li and CI that will eventually collide onto 
the solid-particle surfaces with increasing frequency, and thus could transfer this 
energy {+Q) back into the solid phase to initiate chemical dissociation of LiCl 
solid. These two thermal-assisted processes (due to ±Q) are the causes for this 



16 



ANDREW DAS ARULSAMY 




Figure 3. Both atomic Li and CI in the Hquid phase are denoted 
with fiUed red circles, while the blue circles represent the ionic Li 
and CI in the solid phase. The two blue islands labeled A and 
B are solid particles, immersed in the liquid-phase. This system 
has been deliberately isolated at the critical point when T = Tg 
= 610°C such that the solid phase coexists indefinitely within the 
liquid phase. The magnified sketch beneath the main diagram 
shows the temperature differences between the solid phase (Tgoi) 
and its surrounding liquid phase (Tuq) such that Tuq > Tg > Tgoi 
and Tiiq — Tsoi Tg. At the boundary, the system is in extreme 
nonequilibrium where Li or CI from the liquid phase may react to 
form a rigid ionic bond, and release heat (— Q) as energy. Con- 
versely, the same amount of heat is absorbed {+Q) by Li or CI in 
the solid phase so as to break free from the solid to liquid phase. 



thermal-assisted (or finite-temperature) ^QPT. In a macroscopic point of view, this 
particular isolated system, containing both solid and liquid phases, is in equilibrium 
because the average rate of melting and solidification is the same, hence the solid- 
to-liquid and liquid-to-solid transitions (^QPT) are in balance. But we need to go 
deeper to track the physico-chemical processes at the solid] liquid interface, which 
can be done with lET. For instance, the energy level spacings in the liquid and 
solid phases are Cuquid ^^'^ ^soM' respectively, and since the valence electrons in 
the liquid phase are all in the excited states (thermally polarized) then this implies 
^Mquid < ^soiid- ^his inequality is strictly valid because the excited energy-level 
spacings are always narrower due to weak electron- electron (e-e) interaction in the 
presence of weak electron-nucleus (e-nuc) attraction. Conversely, large energy level 
spacings are inevitable for the energy le vels c lose to the nucleus p]. 

Substituting Cululd < ^soSd into Eq. Q leads to iiqS'+ > soiS+ and iiqS'+ > 
soiS^ , which then allow one to correctly conclude S'nquid > S'soUd where S'uquid — 
liqS'^ + liqS^ and S'soiid = soiS'^ + soI'5'd^- However, as we have noted earlier, 
'S'soiid < 'S'oqin < S'liquid IS invalid because we need to take the nonequilibrium effect 
(at the solidjliquid interface) into account. We now know that this effect occurs 
maximally at the critical point (when T = Tg), and at the solidjliquid interface 
where 5*,"°^ = S'uquid + Ssoiid + S^^lf, and therefore S'soiid < Suquid < S^°^ that 
guarantees the existence of a first-order TPT, while Sf"*g'' on the other hand, ensures 
the existence of a ^QPT, as well as the coexistence of both solid and liquid phases 
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(see Fig. [3]). This means that if S'l^cf 
have S^°^ ^liquid for T > or S","™ 
non-divergent and renormahzed 



— 0, then T ^ Tg and consequently we 
S'soiid for T <Tg. Using Eq. the 



(5.1) C,iTg,t) 



non 
cqm 



a(r)exp<^ 



•^j'^solid + (1 ^ >^j) liquid 



in which, we have defined ^(t) = J2j -^jCsoiid+(l--^j)6iquid where j = {<i,i2, ■ • • ,in}, 
J = X^^i^^/X^^fJ^i and J G [0, 1]. Here, X^^j^^ is the number of Li and CI atoms in the 
sohd phase only, and yi^^?li total number of Li and CI atoms in both liquid 

and solid phases (or in a given system). Moreover, J does not necessarily increase 
with time, if the heat-exchange fluctuates (iQ), and the total time between ti and 
tn is the ti me t aken for a complete solidification of liquid LiCl for constant T — Tg. 
Equation (5.1 ) strictly implies that the magnitude of C't)(t)oqm does not change with 
T only, but also with respect to changes in the interaction-strength parameter, ^(t). 
In particular, Cvit)eqm can change due to changes in ^(i) from ^nquid to ^soUd during 
solidification of LiCl liquid. In this latter case, the heat-exchange as depicted in 
Fig. [3] is solely used to change the interaction strength via ^soiid ^ ^liquid ■ 

If one employs an unrenormalized specific-heat equation, then it is always di- 
vergent for T — Tg because it is undefined at this critical point. On the other 
hand, Eq. (5.1) is well-defined such that ^ can be exploited at the critical points, 
without any divergence. For example, ^liquid and ^soiid are constants for T > Tg 
and T < Tg, respectively, and any heat exchanges that may exist between a given 
system and its surrounding only decrease or increase the system's temperature. At 
the critical point however, the heat exchanges are used only to modify the system's 
physico-chemical properties, hence, the system's temperature remains constant. In 
other words, the renormahzed specific- heat equation (Eq. ( |5.1[ )) reduces to 



(5.2) 



(^liquid^j.) ^ Ci'^^id^T) 



exp 



for liquid LiCl (T > Tg), and for solid LiCl {T < Tg), one just need to switch 
the label "liquid" with "solid" in the equation above. However, neither of these 
equations can be applied when T = Tg. We need Eq. ( |5.1[ ) for T = Tg. This 
completes the proof for Statement 1. 

If we now follow the horizontal arrow given in the inset of Fig.|3j we can show that 
both transitions (for vertical and horizontal arrows) are thermal-assisted ^QPT, 
one occurring during solidification (discussed above for the vertical arrow) and the 
other originates due to changing chemical composition (at a constant temperature). 
For example, the first j^QPT at the critical point during solidification (constant 
T ~Tg) is initiated by the heat-exchange between the liquid and solid LiCl, where 
dS = 5Q/Tg, and from Eq. p^O]), 



(5.3) 



5Q 

Tg 



{Eo ± dC) 



Tg 



15^1 e 



Equation (5.3 1 correctly implies: (a) increasing absorption of heat in a system 



increases the entropy of that system, or vice versa, (b) any amount of change in 
Q (SQ) or ^ (d^) corresponds accordingly to a change in entropy, as it should 
be, however, (c) the entropy of a given system decreases if ^ increases, as strictly 
required by Eq. (3.9), (d) systems with large ^ need large amount of heat to initiate 
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changes such as inching, for example, from Table [2j we know [^Br < Cci] ~^ [^LiBr < 
^Lici] a-nd therefore T^'^'' < T^''~^\ and (e) large amount of heat is required to be 
removed or added to initiate large changes to ^, for example Cuquld ^ ^soUd- Here, 
^Uquid a constant, while Q is the amount of heat removed to initiate the change, 
^Uquid ^ Q ^ isohd- Therefore, (d) refers to the entropy-change of a system prior 
to any phase transition, i.e., for constant ^ and T <Tq oi T > Tg. In contrast, (e) 
reveals the changes in the intrinsic physico-chemical properties of a system due to 
changing ^ when T ^Tg. This means that d\S^\ is only valid at the critical point, 
or when T — Tg because ^ does not change with T, but it changes significantly 
when the physico-chemical properties of a given system change. 

In order to understand the existence of thermal-assisted "-"QPT due to doping 
(follow the horizontal arrow in the inset of Fig.[3|, we increase Br content, replacing 
CI for constant T to obtain a system defined by LiCli_aiBraj . The inequality, 
Cfir < Cci implies ^ decreases with increasing Br content where this doping is carried 
out for constant T ~ 552°C. Now, the critical point can be obtained for ai = 1, 
and at this point, T — Tg = 552°C, which is the melting point of LiBr. As a 
consequence, if —?■ ^^jjfj is achieved through doping at T = Tg ''^\ then Q 

activates the melting process, such that ^^oifj + Q — > ■fuquid- However, we point 



out that Q — <^|^gu[(i — false. Note here that the above stated doping- 



^liquid • 

^liquid ^solid 

induced ^QPT can occur for any constant temperature, and even for T = OK. For 
instance, electronic phase transition (metal-to-insulator) due to doping can occur 
in doped-cuprates at zero Kelvin, giving rise to the usual continuous quantum 



phase transition (okQPT), as explained by Sachdev [23]. Nevertheless, Eq. (4.7) 
remains true for both ^QPT and okQPT, in accordance with the proof given for 
the horizontal arrow of Fig. [ij 

Finally, let us now try to understand what is the implication of having ^-dependent 



^ on Eq. (2.7). In view of Eq. p.U, Eq. (2.7) reads. 



~>c— ion frri .Nnon _ ( ^)(~^^) 



2 5Z "^J^^oUd + (1 ~ Jj)^^quid^ I 



(5.4) X exp <^ — /irA exp 

j 

and this is the electron-ion potential operator at the melting point. In the solid 
phase, 

/r c\ -rVe— ion /rp .Nnon ^ r/e— ion /tsolid\ 

\°-^) "Coulomb V-' 9 ' '-ieqm ~^ ^^Coulomb'.^A J: 

while in the liquid phase, 

(5-6) V^oZm,iTo,tr,l ^ V^co-rn.b(^A''"'')- 

As such, it is clear that there is a finite-temperature quantum phase transition 
from a liquid to a solid phase for a system consisting of (anion-forming) atomic A 
and (cation-forming) atomic B, due to i/^B + ^ab + VS-ZnJi'r''') ^ H^^m + 
^Couionib(^A'''^) as pointed out in the introduction, and in the Further Analyses I. 

In summary, the so-called activation energy in this chemical association should 
be equal to Q, which is required to complete the liquid-to-solid transition, or vice 
versa. In other words, to complete the ^J^'quj^ — >■ C^'ud CsoM ~^ ^Uquid transforma- 
tion, significant changes to the t-dependent many-body wave function are necessary, 
for example, ^'(i)iiqu'id ~^ *(Osond- ^he existence of such a transformation in \E'(t) 
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has been proven within a new quantum adiabatic theorem developed for chemical 
reactions |32j . which can be used to understand why the wave function of unre- 
acted species need to be combined linearly or written in a different form for the 
compounds formed after the chemical associations. Interestingly, Miiller and God- 
dard |41j have also pointed out such a case must exist during chemical associations. 
Thus far, we have exploited the formalism developed for lET such that the only a 
priori information one required to know is the atomic energy-level spacings listed 
in Table [2 



6. Further analyses II 



Here, we apply Eq. (5.1) to a water-to-ice phase transition during freezing. Fig- 



ure [4] sketches the specific heat versus temperature curve for constant pressure 
{Cp{T)) for water (see the dashed line). The dashed line denotes the usual Cp{T) 
curve with a sudden drop (C^'**'"'(r) > Cp™(T)) at the freezing point. Whereas, 
the solid line in Fig. [4] captures the whole mechanism of phase transition from wa- 
ter to ice, including the freezing curve occurring at the freezing point (273. 16K). 
Obviously, the freezing curve is always hidden at the freezing point, and appears as 
a sudden drop in Cp(T) measurements, when the temperature is lowered. Unfor- 
tunately, we do not know the explicit microscopic equation for Cp(T) = dH/dT|p 
where H is known as the enthalpy. However, the fun dame ntal energy-level spacing 



renormalization factor (the exponential term in Eq. (5.1)) for Cv{Tg,t) 
remain intact for CpiTg^t) 



)-- should 



such that 



(6.1) 



Cp{Te,t) 



non 
cqm 



Cp(T)exp. 



Jj^icc ^" (1 



The only change in the renormalizing factor is the numerical constant, 3 in 3 A/2 (see 
the exponential term in Eq. (|5.1[ )). The constant 3 here is unknown and therefore, 
we replace it with rj in Eq. (6.1). This numerical constant, 77 can only be obtained 
by solving Cp{T) = dH/dT|p, but it is irrelevant here. What we actually need are 
the renormalized specific heats for water and ice. 



(6.2) 



3.3) 



CICC 

P 



(T) 



C-'=(r)exp 



A6c 



In Eqs. (6.2) and (6.3|, only C^^*'=^(r) and C'^'^iT) are T-dependent variables, 
and the rest are just 'J -independent variables or constants. This means that 
Eq. (6.2) captures the specific heat for water, whereas, Eq. (6.3 1 is for ice (see 



the dashed and solid lines for water and ice in Fig. |4]). The sketched freezing curve 
however, follows Eq. (6.1), and is smooth here, implying there is no fluctuation in 
Q during freezing, i.e., there is a continuous extraction of heat from the water-ice 
system. Recall here that the freezing curve in Fig.|4]is not observable from the Cp 
versus T measurements alone, but exists as a sudden drop in the specific heat for 
T = 273. 16K (freezing point) because C™'<^''(r) > C^'='=(T) or C""^" < f This 
is the reason why temperature is not a proper variable to monitor at the freezing 
point of any thermally-driven finite-temperature continuous quantum phase transi- 
tion. To those who are not comfortable with this reason, we aggravate the situation 
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Figure 4. A dashed line is sketched to represent the changes of 
constant-pressure specific heat {Cp(T)) with temperature for both 
water and ice. The sohd hne captures the same specific heat based 
on Eqs. (6.1), (6.2 1 and (6.3 1, which gives a complete picture. A 
complete picture here means the ability to capture the changes 
to Cp(T) in the water and ice phases, as well as during freezing 
via £,{t). We have identified temperature T as T{t) because the 
change in temperature corresponds to a change in time. Similarly, 
we also have £,{t) to handle the time-dependent fiuctuations at 
the freezing point (273. 16K). The freezing curve sketched above is 
smooth. The freezing curve will not be smooth if the heat is not 
extracted continuously from the water-ice system, implying the 
heat extraction fiuctuates with time. 



by claiming — in any quantum matter, Eq. (6.1 ) proves that there is no such thing 
as Landau first-order transition. 

In summary, we have proven that ^(t) is the proper quantum variable to capture 
what is really happening at the freezing point. For example, at the freezing point, 
the heat that is being extracted from the water-ice system, activates the energy-level 
spacing transformation, — ^ ^'^^^ or the wave function transformation vj/^^.t^'' — >. 
^icG -j-p ^Yie formation of permanent hydrogen bonds. This transformation 
resuhs in C™^'°''(T) > C^™(r) that implies C*^*"' < f'^^ The latter inequality has 
its origin in hydrogen bonds. For example, both ^'"^.tor ^^^^ ^ico j-gfgj. energy- 
level spacings due to hydrogen bonds, and not due to covalent bonds between O 
and H in a H2O molecule. Obviously, the energy-level spacings due to hydrogen 
bonds are larger for ice compared to water because the hydrogen bonds in ice is 
permanently bonded (static) and therefore, its strength is larger in ice than in 
water. The hydrogen bonds in water phase is dynamic such that the bonds are 
broken and formed randomly. As such, indeed ^'"'^.t®'' < Now for the water 

vapor, all H2O molecules are isolated, and they do not form any hydrogen bond. As 
a result, ^™p°'' refers to the energy-level spacings due to covalent bonds (between 
O and H) within a H2O molecule. Since covalent bonds are much stronger than 
hydrogen bonds, one has ^w^*'''' < ^1^=° < .cvapor ^j^^^^ correctly corroborates with 
the experimental observations [Hiilil], C*''*'=^ (4.187 kJkg^^K-i) > C],™ (2.108 
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kJkg-iR-i) > C™P°'' (1.996 kJkg-iK-i). This concludes our analytical proofs 
and technical analyses on H2O system, again supported by experimental results. 



7. Conclusions 

We first entered into the realm of thermodynamics and statistical mechanics to 
show that all thermal phase transitions occurring in any matter is a proper subset 
of quantum phase transition, even though at times, it is convenient to ignore the 
quantum effects, for example during the transition from a solid to liquid phase. 
Wc went on to rcnormalize the entropy and the Bosc-Einstcin statistics solely to 
prove the existence of the non-divergent specific heat capacity at the critical point 
via the energy level spacing renormalization method. At the critical point, the 
thermal energy is used to change the energy-level spacing of a system, which in 
turn changes the physico-chemical properties of that system. Those properties are 
of course different if the system is in solid phase, compared to liquid. The critical 
points we chose to investigate were the melting points of single-valent ionic bonded 
salts. 

Here, we have formally shown that the finite-temperature continuous quantum 

phase transition (^QPT) is responsible for the thermal phase transition during 
melting and solidification processes, as well as during substitutional doping at a 
constant temperature. Along the way, we also have proven that the first-order 
transition in the constant-pressure specific heat between water and ice satisfies a 
^QPT. In view of the analyses for alkali halide salts and water-ice system, we 
found that the ^QPT is responsible for the liquid^solid and water-ice phase 
transitions, supported by experimental results. In addition, we have given sufficient 
proofs to claim that all first-order thermal phase transitions are equivalent to finite- 
temperature quantum phase transitions, which will lead us to the unification of 
thermodynamics and quantum mechanics. 

The formal proofs developed herein allow one to obtain the precise physical 
mechanisms involved right at the melting point of a given solid, beyond the usual 
physics of a discontinuous thermal phase transition. In other words, our approach 
avoids the diverging thermodynamic variables (when T — >■ T^) at critical points 
(Tc), which are usually in the form of 1/(T — T^Y where o is some constant. For 
example, wc have been able to formally prove the existence of a finite-temperature 
quantum phase transition due to some wave function transformations (giving rise 
to changing ionization energies) during the melting of alkali halides, and also when 
water freezes over. This means that one can gain a far more detailed understanding 
of the physical phenomena occurring at a critical point as a result of quantum 
effects (namely, wave function transformations that give rise to a quantum phase 
transition). 

In particular, we have shown that the discontinuous thermal phase transition is 
actually a continuous quantum phase transition, which is observable as an electronic 
phase transition. This electronic phase transition is the cause for the solid-to-liquid 
(or liquid-to-solid) and liquid-to-gas (or gas-to-liquid) transformations with well 
defined (non-divergent) thermodynamical properties at the critical points. More- 
over, wc found that the theoretical predictions extracted from the proofs are very 
well supported by the experimental results obtained from the alkali halides, namely, 
the atomic ionization energies, melting points and diatomic bonding energies. For 
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water and ice, we had the constant-pressure specific heats and atomic ionization 
energies in complete agreement with the time-dependent specific heat formula. 
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8. Appendix 



1. Renormalization of Eq. (2.4). The renormalization of Eq. fl2.4| is straight 



forward where we substitute all k (interaction potential constant (also known as 
the spring constant)) with Eq. (2.8) to obtain 



(8.1) 



2 



'fc- (e727reoi?3) jk + {e^ /2TreoR^) 



m \j m 

Any variable found to wear a tilde, here or in the text, is a renormalized parameter. 



8.2. Derivation of Eq. (2.6). We start from Eq. (8.1) and expand it in series in 
the form of ^Jl + x = l + x/2- + x'Vl6 , 



E = 

2 



'fc e2/27reoi?3 

TO TO 



TO TO 



(8.2) 



1^ /A 
2 V TO 



e2/27reoi?3^ 



k 



e2/27reoi?3 



Substituting Cjq = \lk/m into Eq. (8.2) and noting that x — ^{e^ /2TreoR^k), 

2 



E — -hujn 
2 ° 



(8.3) 



1 



1 



^^^[^ 27reoi?3A:) 8( 27reoi?3A:) ^ 2Tre„R3k 



1 



2V27reoi?3fcy 8\2TT€oR^k 



16 V27reoi?3fc 



It is clear from Eq. (8.3) that all the terms with odd powers {x^-^'^'"') cancel out, 
and the terms with only even powers, x^'^'^' " survive. Consequently, we just need to 
impose the condition, A;exp[A^] [e^ /2'KeQR^), which means {e^ /2T:eQR^Yk'^ « 
for all n > 4 where n € Neven- Here, Ncvon is the set of even natural numbers. After 
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doing SO, 



(8.4) 



E = -hiUn 

2 ° 



1 - 



1 - 



To see if there is any (repulsive or attractive) interaction between the two identical 
and neutral atoms (see Fig. [2](A)) as a result of Hi, we need to subtract hujQ (the 
solution for Hq) from E (the solution for Hq + Hi) such that 



(8.5) 



E{^) — hujQ — hujQ 
= hujQ 



2- 



1 / 



4fc2 V27reoi?3 
P \2TTeoRy 



- 1 



We now substitute = loqW? and then followed by the substitution 9J, tZiQ 
woexp [^A^] into Eq. (8.5) to arrive at 



,6) E{() - m = VSi,.(0 = 



exp 



The negative sign in Eq. (8.6 1 implies there is an attractive interaction, and this is 



the van der Waals attraction in its renormalized form. But Eq. (8.6 1 is only valid 
for small polarization, made sure by the condition [37], A:exp[A^] » (e^/27reoi?'^). 
For large polarization, which is necessary at melting points, we take fcexp[A^] w 
(e2/27reoi?^) and this leads us to rewrite [37] Eq. (fsjj). 



(8.7) 



E{0' = 2^ 



Hence, the van der Waals attraction in the presence of large polarization [37^, 



1 / Oh 



Again, noting that k — fcexp [A^] and Cjq — ujq exp [^A^], we can simplify Eq. (8.8), 

(8.9) E{0' -fuLo = fi^o '' 

Since the attractive interaction is mutual between two identical and neutral atoms. 



we can then average their contributions by rewriting Eq. (8.9), 
(8.10) ^(O'-^^o = 













j 1 exp 




+ exp 


^A^B 


} 



and this is Eq. (2.6). 
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8.3. Derivation of Eq. (2.5|. The attractive interaction captured in Eq. (8.10) 



is due to large polarization in the absence of Coulomb repulsion between the po- 
larized electrons |37| . If we inclu de the Coul omb repulsion, then we strictly require 
V^Couiomb < |V^waais(C)l (see Eqs. (Elll and (Ho)) where, 



ioulomb 



.11) 



1 (-e)(-e) 

2 47reo|rA - 

+ exp ( - /i(rA 



exp I 



rB)e 



Now to guarantee an attractive interaction without the additional condi tion, Vcouiomb < 
l^^aais(C)l: uccd to eliminate the Coulomb repulsion given in Eq. (8.11) by im- 
posing an asymmetric polarization as depicted in Fig.[2]jB). This asymmetric polar- 
ization is obtainable between a cation and an anion where the attraction is between 
the anion (ion A, charge -I-) and the outer electron (e, charge — e) of a cation (ion 
B) [37] such that 



T 7-e— ion 
Coulomb 



.12) 



l[ (-e)(+e) 
2 1 47reo 

, i-e){+e) 

47reo|RA - tb 



1 



|Ra " tb 

1 

R.B - TaI 

exp ( - 



|a5a 



exp [~ fJ-r 
exp ( — ^ree^^Afo^ 



fir AC 



Obviously, the second term in the curly bracket on the right-hand side of Eq. ( 8.12 ) 
is zero and the factor, 1/2 disappears because the attraction now is unidirectional. 
This means that the outer electron of the easily polarizable cation (B) is attracted 
to an anion (A), which has the least polarizable electron, i.e., the Coulomb repulsion 
has been eliminated. Putting Eqs 



101 and (8.12) together gives Eq. (2.5 1 



8.4. Derivation of Eq. (3.8). Technically, it is easier to work with sums 
compared to products (J^^JTand therefore, we convert 



.13) niN,q)=l[-j^^^^^lnniN,q) = Y^ln- 



ni\{qi-n^)l ' ^ ni!((ji - n^)! ' 

As a consequence, 

(8.14) lnr!(7V,g) =^lng,!-lnn,!-ln(g,-n,)!. 

i 

Note here that this conversion does not disturb the "exactness" of the derivation in 
any way because we are interested in the maximum of fl{N, q), not in its absolute 
value such that 

(8.15) d[\nn{N,q)] ^^^ dm.q) ^ 

dui drii 

Using Stirling's approximation, we can write InnJ = rii Inn^ — n^, Inqil — 
J2i Qi In Qi - Qi and ln(gj - "»)! = J^Mt - ^i) ln(<?i - Ui) ~ (g^ - n^), substituting 
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these into Eq. ( |8.14[ ) will lead us to 

lnf^(A^, q) = ^ {qi Ing., - (?i - (n^ Inn., - n^) - [{q, - iii) \n{q, - Ui) - {q^ - n,;)]} 

i 

(8.16) = {qi In qi - rii \n Ui - qi\n{qi -Hi) +ni\n{qi-ni)}, 



The next step is 

d[\nn{N,q)] 
drii 



1 


h In Hi 




-1 






n, 




qi 


f 


+ 


Tli 






qi-rii 







E 

^ < - 1 - lnn.i + \n(qi 



qi - n. 



+ ln(gi - Hi) 



qi - ni qi - Hi 



i.l7) - ^In 



qi - n. 



Hi 



We now substitute Eq. ( |8.17| ) into Eq. 

dFncw sr^, qi — ni ^^-^ dui v^/cn , /-\ 



drij 
dui 



(8.18) 



qi - rij 



After summing over all the particles (rij) and quantum states {qi), one gets 



In 



(8.19) 



q- N 
N 
N 



N 



q~ N 

cxp [n + \b.{Eq ± ^)]. 



exp [^1 + \b{Eo ± Cj], 



N{qlN-l) 

Equation ( |8.19| ) leads straight to Eq. (3.8) after replacing /i = — /i and Ab = — Ab. 
These replacements are necessary to make sure lim^_s.oo limAg_i.oo N /q ~ 0. 

8.5. Derivation of Eq. ( 3.16[ ). The derivation here is similar to Eq. (3.8), hence 
we can skip the conversion (J^^ — ^ In) and the Stirling's approximation, and 
start from 

5[lnD(7VD,gD)] 



d 



drii 



^{rii ln(ni + qi) + q^ \n{ni + qi) - Ui Inn,; - qi IngJ 



E 



+ ln(ni + qi) 



rii + qi 

i ^ 

(8.20)= El-^- 



1 



'ni + qi 

qi 



rii h In Tli 

n, 



rii + qi ni + qi 



We now recall Eq. ( |3.15p 



.21) 



dm 



E\ , n-i + q., dui r ly^.. 



drii 
dui 



= 0. 
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After summing over, 
(8.22) In^^+a 



Subsequently, we get 
(8.23) 



^\B[Ee^^^] = 0. 
exp [a + Ae^^e 5 



Equation (8.23) is nothing but Eq. (3.161, after replacing a — —a and Ab 
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